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We show the spin-orbital coupling induced by an artificial hght-induced gauge field can fully 
restore superfluidity suppressed by population imbalance in a two-dimensional (2D) Fermi gas, 
leading to unconventional superfluid states either with topological Majorana fermion excitations 
or showing a novel mixture of triplet pairing with spin-up (down) components respectively in the 

± ipy pairing channels. We self-consistently calculate the zero temperature phase diagram at the 
BCS side of Feshbach resonance and show that the phase transitions between different superfluid 
states can be revealed through measurement of the in-situ density profile of the 2D atomic cloud in 
a weak global trap. 
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Ultracold atomic gas provides an ideal platform to 
study superfluid states under tunable configurations. 
There has been great interest in looking for unconven- 
tional superfluid states beyond the conventional s-wave 
superfluidity observed in experiments [H . Population im- 
balance provides a mechanism to suppress the s-wave su- 
perfluid state and can lead to novel superfluid phase, such 
as the FFLO state . Polarized Fermi gas with tunable 
population imbalance has been studied extensively both 
experimentally and theoretically [Hill- The FFLO state, 
unfortunately, is typically fragile and hard to observe 
experimentally P, Increasing population imbalance 
usually destroys the superfluid state instead of leading 
to novel superfluidity. 

The spin-orbital (SO) coupling induced by an artifi- 
cial gauge field emerges as a powerful control method for 
ultracold atomic gas d-Q. Various proposals have been 
made to realize SO coupling in ultracold atomic gas 
and remarkable experimental progress has been reported 
to demonstrate the artificial gauge field @ . Motivated by 
this progress, strong interest arises recently in studying 
the three-dimensional ultracold Fermi gas under artificial 
SO coupling [6j]. It has been noted that the SO coupling 
can enhance superfluidity for this system [6] . 

In this paper, we show that interesting and exotic 
superfluid states arise in two-dimensional (2D) Fermi 
gas from interplay of population imbalance and SO cou- 
pling. Without SO coupling, population imbalance above 
a critical value suppresses superfluidity, leading to nor- 
mal phases. We show that SO coupling fully restores 
superfluidity in 2D, but the pairing is not in the con- 
ventional s-wave channel any more and the resultant su- 
perfluid phases have more exotic properties. Under large 
polarization, one of the spin components dominates with 
the pairing in the + ipy channel. Such a phase sup- 
ports non-abelian topological excitations by trapping a 
Majorana Fermion in its vortex core [7|. Thisphase has 
been predicted before for a 2D Fermi gas 0, [1] , and the 
contribution of our calculation of a self-consistent phase 
diagram is to determine the stability region of such a 



phase under real physical parameters, which is impor- 
tant for experimental observation. Furthermore, our cal- 
culation predicts a new kind of unconventional super- 
fiuid phase under intermediate polarization, with spin- 
dependent pairing in the p^iipy channels. In this phase, 
the majority (minority) components pair up respectively 
in the Px + ipy {px — iPy) channels. We characterize the 
transition order between different superfluid states, and 
show through explicit calculation that these transitions 
can be revealed by observation of the singularity points of 
the density profile of a 2D atomic cloud in a weak global 
trap, which can be measured in experiments through the 
in-situ imaging Q. 

In experiments, one realizes the 2D Fermi gas by ap- 
plying a strong optical trapping potential (or an optical 
lattice) to ultracold atoms along the ^-direction. Near 
the potential minimum where the atoms are located, the 
potential is well approximated by a harmonic trap with 
Vz {z) = 7710^2^:^/2, where m denotes the atomic mass and 
LOz is the trapping frequency. At the BCS side of the 
Feshbach resonance, the transverse wave function for the 
atoms is given by the ground state of Vz (z), and the 
atomic collision can be characterized by an effective 2D 
interaction [iG]. With light induced Rashba-type of SO 
coupling [^-[a], the Hamiltonian for the system can be 
described as 

H =^ Ck,<Ta[,CTak,o-+ k |^e"^'''a|^ ,f.ak,; + h.c. 

k,cr k 

+ ^ X! 4,t«-k-Hq,ia-k'+q,iak',t- (1) 

k,k',q 

where Ok.o- and ^ denote the fermionic field opera- 
tors with in-plane wave vector k = {kx,ky) and spin 
a =ti 4" The free particle dispersion relation is given by 
Ck,a = Ck - M(T, where ek = ?i^k^/ {2m) and ^tJ- ^ fJ-^h, 
with fi being the chemical potential and h being the effec- 
tive Zeeman field (equivalent to a population imbalance 
between spin-up and spin-down components). The sec- 
ond term of H describes the Rashba-type SO coupling. 
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whose strength is denoted by the coefBcient a, with k and 
ipk being the magnitude and azimuthal angle of the the 
in-plane wave vector k. A combination of the Zeeman 
field and the Rashba-type SO coupling can be realized 
through control of a few laser beams |^]. The Ub term in 
the Hamiltonian describes the effective 2D interaction, 
where L is the quantization length and Ub is the bare 
coupling rate. The bare coupling Ub is connected with 



the physical coupling 
scattering length [lOl 
relation 



Up {Up is determined by the 3D 
) through the 2D renormalization 
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(2ek + Jii^z 
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We determine the self-consistent phase diagram of 
the Hamiltonian H within the mean-field framework, 
which is a reasonable approximation at zero tempera- 
ture. Under this framework, we introduce the pairing 
order parameter A = C^f)i~^ X^k (^-k,iak,t) to decom- 
pose the last term of H into '^t -f^-k i + H.c^ — 

\A\^ L'^ /Ub through the Wick theorem. The quadratic 
mean-field Hamiltonian can then be diagonalized into 

H^Y. (^k,+4,+"k,+ + i;k,-aL-"k.-) +L^^o, (3) 

k 

where Q;k,± denote the quasi-particle modes with the ex- 
citation energies 



+a^k^ + h^ + \Af±^J{h^ + a^k'^)^l + h'^\A\' 



,(4) 



and = Ck — M- The last term in Eq. ([3]) denotes 
the zero-temperature thermodynamic potential, with fio 
given by 

^0 = L-^ - (^k,+ + i?k,-) m - \Af /Ub. (5) 



At zero temperature, the quasi-particles dk,± are in vac- 
uum states and we should minimize VIq with respect to 
the order parameter A to determine the ground state of 
the system. 

Usually one calculates the phase diagram of the system 
by solving the gap equation dClo / dA ~ and the number 
equations 



(6) 



However, this is not a reliable approach here since the 
gap equation has multiple solutions corresponding to un- 
stable and metastable states. A typical thermo-potential 
r^o is shown in Fig. [1] As one can see, it has multiple 
minima and maxima satisfying d^^/dA = 0. To find the 
ground state of the system, in this paper we directly cal- 
culate rig for various A to figure out its global minimum. 
In our calculation, we take Tilo^ as the energy unit and 
= \/'h/ (mujz) as the length unit. The units for Up 
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FIG. 1: The thermo-potential Slo as functions of the su- 
perfluid order parameter A under different effective fields h. 
Other parameters include the chemical potential /i = 0.5, the 
spin-orbital couphng strength a = 0.1, and the interaction 
parameter Up = —5 (corresponding to a 3D s-wave scattering 
length as ~ ~2az). Normalized with the energy unit fiu)z and 
the length unit a^, all the parameters are dimensionless in 
this and the following figures. The arrows indicate the global 
minimum of the corresponding curves. 
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FIG. 2; The superfluid order parameter A (a), the number 
density n (b), and the magnetization = (n^ — n^)/2 (c) as 
functions of the effective Zeeman field ft under different spin- 
orbital coupling strengthes. The chemical potential and the 
interaction parameter are fixed at = 0.2 and Up = —5. The 
inset of Fig. (a) shows the exponential decay of the superfluid 
order parameter with increase of the Zeeman field. 



and OL are given respectively by a^hLOz and ayhiUz- With 
these units, all the parameters become dimensionless. 

In Fig. [2l we plot the superfluid order parameter A, 
the number density n, and the magnetization Sz = 
{n^ — n^) /2, as functions of the effective Zeeman field h 
under various values of the SO coupling rate a. At a = 0, 
it is clear there are two phase transitions induced by in- 
creasing the field h. Below a critical field he ^ 0.1, we 
have a superfluid phase with zero magnetization. This 
corresponds to the conventional BCS state. Above he, 
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FIG. 3; The phase diagram of a polarized 2D Fermi gas un- 
der SO coupling, with the parameters f/p = —5 and a = 0.1. 
The SF, TSF, SFp±ip, VAC, represent, respecivtely, the con- 
ventional superfluid, the topological superfiuid, the spin- 
dependent mixture of p^, ± i-py superlfuid, and the vacuum 
phases. When the order parameter A gets extremely small, it 
is very likely that the superfluid will be destroyed by thermal 
or quantum fluctuation, leading to normal mixture (MN) or 
polarized normal (PN) phase. The dashed curve represents 
a crossover line by setting A = 10~^ (the choice of 10~^ is 
arbitrary, however, the line has no significant change by set- 
ting a different small cutoff value for A due to its exponential 
decay) . 
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FIG. 4: The occupation number rih (a-c) and pairing func- 
tions (o_k,aak,CT') (d-f) as functions of for different phases: 
(a,d) the conventional superfluid phase with parameters [i — 
0.3, h — 0.1; (b,e) the spin-dependent mixture of px ± ipy 
superfluid phase with ^ = 0.3, h = 0.2; (c,f) the topological 
superfluid phase with /i = 0.05, h = 0.1. 



the superfluid order parameter A is suppressed to zero, 
and we have a finite magnetization Sz with rif > > 0. 
This corresponds to the normal mixture phase. Further 
increasing the field above he ^ 0.2, the minority compo- 
nent nj. is reduce to zero and we have a polarized normal 
phase. 

With a finite a, the picture changes qualitatively. At a 
small a, such as a = 0.1 in Fig. 2, there is still a big drop 
of the the superfluid order parameter A when the field h 
is above a critical value h^, but A does not drop to zero. 
The critical field he increase with a, which indicates that 
the SO coupling enhances superfluidity. The transition is 



still of the first order for small a similar to the a = case. 
What happens is that the thermo-potential VIq as a func- 
tion of A has two non-zero minima with A2 > Ai > 0, 
as shown in Fig. 1. As h increases, Ai, replacing A2, 
becomes the global minimum, and the ground state un- 
dergoes a first-order transition between different types of 
superfluid phases when the order parameter jumps from 
A2 to Ai. Above the critical field he, the order param- 
eter A eventually decreases exponentially with increase 
of the field h as shown by the insert of Fig. 2(a), how- 
ever, it does not reach exact zero within the mean-field 
framework. Due to the exponential decrease of A, the 
order parameter becomes very small for large h and it 
will be destroyed by thermal fluctuation even at very low 
temperature or likely by quantum fluctuation beyond the 
mean-field framework. At larger a, the sudden drop of 
the order parameter A disappears, as indicated by the 
a = 0.2 case in Fig. 2(a). If we look at the derivative of 
A as a function of ft, at a = 0.2, it shows a kink at the 
critical he, suggesting that the phase transition changes 
from the first order to the second order for large a. 

With nonzero SO coupling, the superfluid phase above 
the critical he has novel features and is the focus of our 
interest in the following discussion. We concentrate on 
the small a case, as it is easier to realize this case in ex- 
periments. Below the critical he, the superfluid phase is 
similar to the conventional BCS state, where the singlet- 
pairing dominates. Above the critical he, there are two 
kinds of unconventional superfiuid phases. From the or- 
der parameter A, it is hard to distinguish these phases 
as A changes continuously across the phase boundary. 
However, from the number density and the magnetiza- 
tion shown in Fig. 2(b) and 2(c), one can clearly see a 
kink at the boundary, indicating a second order phase 
transition. The boundary is determined by the condition 
/i^ = /i'^ + A'^. In different superfiuid phases, the excita- 
tion spectra i^k,± of the quasiparticles are always gapped, 
except at the critical point with h = \/ -\- A^, where 
the excitation spectrum -Ek,- becomes gapless. When 
h > yTi^ + A^, the phase is identified with the topolog- 
ical superfluid state which supports Majorana fermion 
excitations with exotic n on-abelian fractional statistics 
0, g. When he< h < y^JJTA^, we have a new kind 
of unconventional superfluid phase whose nature will be 
studied below. 

The overall phase diagram is shown in Fig. 3 at 
a = 0.1. To understand the nature of different super- 
fluid phases, we show in Fig. 4 the corresponding occu- 
pation number (nk,o-) in the momentum space and the 
pairing functions (a-k,o-flk,cr')- In the conventional su- 
perfluid phase denoted by SF there, the pairing is dom- 
inantly in the s-wave channel between the spin t and 
J, components with small hybridization from the triplet 
pairing (a_k,tak,t) and (a~k,4.ak,i) induced by the SO 
coupling. The momentum distribution (nk.o-) is rounded 
off by the strong pairing interaction and (j^k.t) ^^^d (n.k,^) 
are almost identical as the magnetization is small. In the 
topological superfluid phase denoted by TSF there, the 
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FIG. 5: The density profile (left column) and the super- 
fluid order parameter (right column) as functions of the ra- 
dius r± of the atomic cloud under different population imbal- 
ance: (a,c) h = 0.1 and the chemical potential at the trap 
center no — 0.20, corresponding to the total atom number 
N = lO'' and the polarization P = (Nf - N^)/N ~ 0.31. 
(b,d) fio = 0.19 and h = 0.15, corresponding to = 10* and 
P = 0.93. The other parameters are Up = —5, a — 0.1, and 
the aspect ratio A — 400. 



majority spin component (f^k,t) dominates. The pair- 
ing is mainly in the triplet channel (d_k,-|-ak,t) with a 
phase factor of e^'^'' , where (pk denotes the azimuthal an- 
gle of k. This corresponds to the Px + ipy superfluid 
phase, and it is well known that this phase has Majorana 
Fermion excitations with non-abelian statistics [3,S]- For 
the superfluid phase with he < h < y^i^ + A^, we have 
significant population in both spin components. The 
pairing is dominantly in the triplet channel (a-k,t^k.t) 
and (a_k,4,ak,4,)- However, the symmetry is different for 
the tii components: the pairing phase for (a-k,tak,t) 
((a-k.^ak,;)) is given by e"^'= {e-"^"), suggesting that 
the spin-up and down components are in different p-wave 
superfluid states, with px + ipy {px — iPy) symmetries, 
respectively. The system is in a mixture of these two 
spin-dependent p ± ip states, and the phase is denoted 
by SFp±jp there. In a vortex of this superfluid phase. 



both components may support Majorana fermions, but 
a combination of them gives topologically trivial excita- 
tions. As a mixture of two topological superfluids, we ex- 
pect the SFpiip phase has topologically non-trivial prop- 
erties, which could be activated through spin-dependent 
manipulation and detection of the atomic cloud [4j . This 
deserves further study in future. 

To have experimental signature of the transition be- 
tween different superfluid phases, we look at the den- 
sity profile of the atomic gas in a weak global harmonic 
trap, which can be measured directly in experiments Q. 
Assume the global trap in the x-y plane is given by 



V± {x, y) = vnjjj\ (a;2 y'^) /2 



/2 with the as- 



pect ratio A = ijJz/ijj±^ 3> 1. The density profile can be 
calculated through the standard local density approxi- 
mation which replaces /i in the homogeneous case by the 
position dependent ^i{r±_) = — V_\_(x,y), where 
denotes the chemical potential at the trap center. The 
typical density profiles and the corresponding superfluid 
order parameters are shown in Fig. 5 under two scenar- 
ios. In Fig. 5(a) and 5(b), we see a first-order transition 
from the conventional SF phase to the topological TSF 
phase, signified by a jump in the density profile for both 
(nk.t) and (rik,i}- In Fig. 5(c) and 5(d), we see a second- 
order phase transition from the SFp-t^p phase to the TSF 
phase, indicated by a kink in the density profile of (nk,i)- 

In summary, we have self-consistently calculated the 
phase diagram of the 2D polarized Fermi gas under SO 
coupling, and found two kinds of unconventional super- 
fluid phases with either topological non-abelian excita- 
tions or spin-dependent mixing oip±ip superfluid states. 
The transition order between different superfluid phases 
is specified and the in-situ measurement of the density 
profile of the atomic gas provides a convenient method 
to reveal the phase transition and its order. 
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